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ABSTRACT 



> 
CO 

O We analyze the spectral properties of driven, supersonic compressible magne- 

tohydrodynamic (MHD) turbulence obtained via high-resolution numerical ex- 
periments, for application to understanding the dynamics of giant molecular 
clouds. Via angle-averaged power spectra, we characterize the transfer of en- 
ergy from the intermediate, driving scales down to smaller dissipative scales, 
Q-f and also present evidence for inverse cascades that achieve modal-equipartition 

levels on larger spatial scales. Investigating compressive versus shear modes sep- 
arately, we evaluate their relative total power, and find that as the magnetic field 
strength decreases, (1) the shear fraction of the total kinetic power decreases, and 
(2) slopes of power-law fits over the inertial range steepen. To relate to previous 
work on incompressible MHD turbulence, we present qualitative and quantitative 
measures of the scale-dependent spectral anisotropy arising from the shear-Alfven 
cascade, and show how these vary with changing mean magnetic field strength. 
Finally, we propose a method for using anisotropy in velocity centroid maps as 
a diagnostic of the mean magnetic field strength in observed cloud cores. 



1. Introduction 

The molecular interstellar medium (ISM) is observed to be turbulent at all spatial scales 
- from the dense, small prestellar cores where the turbulent Mach number Ai s = o~ v /c s is 
order-unity, to the low (mean) density, giant molecular cloud complexes (GMCs) where M. s 
is 10 or more (e.g. Blitz (1993)); here o~ v is the velocity dispersion and c s is the sound speed. 
Such high Mach numbers imply a strongly compressible state on large spatial scales, which 
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is manifest in the inhomogeneous structure of observed GMCs. Although magnetic field 
strengths are observationally difficult to obtain (see e.g. Crutcher (1999)), magnetizations 
at up to kinetic-equipartition levels are possible for GMCs that form in a flux- conservative 
way from the diffuse ISM. GMCs may thus be thought of as gaseous entities pervaded by 
compressible magnetohydrodynamic (MHD) turbulence. This turbulence is believed to have 
many important effects on cloud structure and evolution, but remains only incompletely 
understood. 

In recent years, large-scale numerical MHD simulations have been introduced as a way 
to explore the fundamental properties of compressible MHD turbulence and hence interpret 
the dynamics of GMCs (see e.g. Vazquez-Semandeni et al. (2000)). An initial emphasis 
of these studies was to assess the dissipation rate of compressible MHD turbulence (e.g. 
Stone et al. (1998), hereafter Paper II; Mac Low et al. (1998); Padoan & Nordlund (1999)), 
with the general conclusion that magnetic fields do not appreciably lengthen the dissipation 
time compared to the ratio L/a v associated in an unmagnetized medium with either the 
crossing time to form shocks (for compressible flow) or the turnover time of an eddy (for 
incompressible flow). Other studies have focused on characterizing various aspects of the 
density, velocity, and magnetic field structure, both intrinsic and as observable with spectral 
or continuum diagnostics (see e.g. Ostriker (2002) and references therein). 

An important property of compressible MHD turbulence that has not previously been 
extensively analyzed is the power spectrum - namely, how the fluctuations of dynamical 
variables (especially velocity and the magnetic field) vary with spatial scale and angular 
direction. 1 These spectral properties reflect how energy is transferred from input scales 
to smaller (dissipation) scales, and potentially also to larger scales. Spatial variations in 
heating may depend on the nature of the "forward" energy cascade (to small scales), while 
the viability of tapping internal sources of mechanical energy (such as stellar outflows) to 
provide large-scale turbulent support of a cloud depends on the efficiency of the "inverse" 
energy cascade. Determining what signatures the input scale leaves on the power spectrum 
may help identify the sources of turbulence. Because spectral properties depend on the 
mean magnetization of the medium, spectral characteristics inferred from molecular line data 
cubes or from polarized extinction/emission maps may be useful for diagnosing the magnetic 
field strength. One of the most well-know characteristics of turbulence in molecular clouds 
is the general power-law increase of linewidth with spatial scale, as initially described by 
Larson (1981). Analytical methods are currently being developed (e.g. Heyer & Schloerb 
(1997); Rosolowsky et al. (1999); Ossenkopf & Mac Low (2002); Lazarian & Pogosyan (2000); 



1 The work of Cho & Lazarian (2002), contemporaneous with our study, also analyzes some spectral 
aspects of compressible MHD turbulence. 
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Ostriker et al. (2001)) that may enable relatively detailed measurements of the velocity power 
spectrum from molecular line observations. 

With the above motivations in mind, in this paper we analyze a high-resolution (up to 
512 3 ) set of supersonic (M s = 5) driven-turbulence, ideal MHD simulations to characterize 
their basic spectral properties. Using four different models with Alfven speeds in the mean 
magnetic field from zero up to 2a v , we focus in particular on how the spectral properties 
depend on the mean magnetization of the system. 

The plan of this paper is as follows: We begin in §2 with a brief review of theoretical 
models and numerical results on power spectra in other regimes of turbulence (incompressible 
and/or unmagnetized), to place the current work in context. In §3 we describe the numerical 
methods used to construct our models and the parameter sets we have adopted. Section 4 
presents the results of our spectral analyses. First, in §4.1, we examine the angle-averaged 
(one-dimensional) power spectra of the magnetic and velocity fields, for shear and com- 
pressive components separately and in combination. Section 4.2 then analyzes the global, 
scale-dependent, anisotropy introduced to the power spectral density (PSD) in the presence 
of a dynamically-strong mean magnetic field. Finally, in §5, we conclude with a discussion 
of our results and comparison to the results from other numerical simulations, theoretical 
model predictions, and current observational measures of turbulence. 



2. Spectral Properties of Turbulence 

Most commonly, "turbulence" refers to a complex state of randomly fluctuating velocity 
fields. Depending on the nature of the medium and the strength of the turbulence, there 
may also be accompanying fluctuations in the local magnetic field and in the density and 
temperature. The scale-dependent structure of the turbulence can be described in many 
ways, with perhaps the simplest characterization in terms of power spectra - i.e. the Fourier 
transforms of the fluctuating part of the fluid variables. In general, a power spectrum -P(k) 
may depend on all three components of the wave vector k. If the turbulence is isotropic, 
then P depends only on the magnitude k = {k 2 x + k 2 . + k 2 ) 1 / 2 of k (also written as k r ), and 
a one-dimensional spectrum averaged over polar and azimuthal angles fully characterizes P. 
When there is a preferred direction in the medium, as for the case with a strong, large- 
scale magnetic field, then turbulence may be anisotropic. For a single preferred direction 
x, P depends only on the components of k parallel and perpendicular to x, and may be 
characterized by a two-dimensional spectrum P(k\\, k±) averaged over cylinders centered on 
the x axis. 
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An equivalent way of characterizing the scale dependence of turbulence is in terms of 
the variance as a function of physical scale; "Larson's-Law" empirical linewidth-size relations 
use this representation. For the case of isotropic turbulence with a power-law spectrum, the 
variance of the corresponding fluid variable will also increase as a power law as a function 
of the spatial scale over which it is averaged. For example, if t> 2 (k) oc k~ n , then the velocity 
dispersion increases as a function of physical size £ as a v (£) oc i^-^l 2 . Larson (1981) 
originally reported a linewidth-size relation a v (£) oc £ - 38 ; larger exponents ~ 0.5 — 0.6 have 
been reported for larger scales with more homogeneous data sets (e.g. Solomon et al. (1987); 
Brunt & Heyer (2002)). Because a number of issues remain unresolved in extracting the 
velocity power spectrum from molecular-line data cubes, at present we may regard observed 
turbulence in GMCs as broadly consistent with a range of possible spectral indices. 

The most familiar model of turbulence is that of Kolmogorov (1941), conceived to 
describe the energy cascade from large to small (dissipation) scales in incompressible, un- 
magnetized flows. Local, isotropic interactions of eddies with velocities v(£) at physical 
scales £ have durations £/v(£). An assumption of a conservative energy cascade within an 
"inertial range" then implies v 3 (£)/£ is scale-independent, so that the velocity dispersion on 
the scale £ is v (£) oc £ 1 ^ 3 . Assuming isotropy and converting to a power spectrum using 
J k d 3 k'P(k') oc v 2 (£) for k = 2n/£, the implied angle-averaged Kolmogorov power spectrum 
thus obeys the scaling P(k) oc /c~ 2 / 3 x k~ 3 oc A;" 11 / 3 . Numerical simulations in the mildly 
compressive limit (decay models with M. < 1, and driven models with M. ~ 1) support 
the velocity scaling predicted in the Kolmogorov model while showing that the compressive 
component of energy amounts to 5-10% of the total (Porter et al. 1994, 1998, 1999). Because 
the Kolmogorov model assumes an incompressible medium, it cannot be expected to apply 
at large scales in a molecular cloud, where the turbulence is strongly supersonic. 

At an opposite extreme from Kolmogorov's energy-conservative, incompressible turbu- 
lence cascade is the Burgers (1974) model of shock-dominated turbulence. In this and related 
models, the presence of shocks transfers energy immediately from large-scale motions to dis- 
sipative scales, rather than via a cascade through "eddies" at intermediate spatial scales. 
The spectrum P(k) oc k~ 4 associated with Burgers turbulence corresponds to the Fourier 
transform of a collection of step functions representing many individual shock profiles. Con- 
version from /c-space to the spatial domain yields v {£) oc £ 1 ^ 2 . 

Neither the Kolmogorov nor the Burgers model incorporates magnetic fields, which 
may be important to the mode of energy transfer through scales - particularly if the total 
magnetic energy density nears that in the turbulent velocity field. Approaches to theoretical 
modeling of magnetized turbulence have focused on the incompressible case, in the limit in 
which fluctuations in the magnetic field are small compared to its mean value; an excellent 
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current review is contained in Cho et al. (2002b). The physical picture in the early theory 
of Iroshnikov (1963) and Kraichnan (1965) (collectively, "IK") is based on energy transfer 
to smaller scales initiated by collisions of counter-propagating Alfven wave packets. An 
assumption of an energy-conservative cascade and isotropy (with the implication that many 
collisions are needed to dissolve a wavepacket) yields a spectrum P(k) oc k~ 7 / 2 . However, 
consideration of the resonant conditions for three-wave interactions among Alfven waves 
suggests that the cascade will in fact not be isotropic, but instead should preferentially 
transfer energy in the direction perpendicular to the mean magnetic field (Shebalin et al. 
1983) (hereafter SMM83). Numerical simulations in 2D and 3D by SMM83 and Oughton 
et al. (1994) directly demonstrated that spectral anisotropy extended in the kj_ direction is 
indeed present; this anisotropy increases toward smaller spatial scales. 

A model of incompressible MHD turbulence that uses the colliding-wavepacket picture of 
IK, while replacing the unphysical assumption of isotropy with the primarily perpendicular 
cascade of SMM83, was proposed by Goldreich & Sridhar (1995, GS; see also Goldreich 
& Sridhar (1997)). GS introduce the concept of a "critically-balanced" cascade in which 
the nonlinear ("eddy turnover") time of a wavepacket and its propagation time remain 
comparable, so that k\\VA ~ k±v(k±). Combining this relation with an energy-conservative 
near-perpendicular cascade (k±v 3 (k±) ~ constant) yields a scale-dependent anisotropy law 
fey oc k 2 / 3 . Numerical simulations in the incompressible limit by Cho & Vishniac (2000), 
Maron & Goldreich (2001), and Cho et al. (2002a) have verified (via structure functions for 
coordinates aligned with the local B-field) that the expected scale-dependent anisotropy is 
present. The power spectrum from these incompressible MHD simulations is similar to or 
slightly steeper than the Kolmogorov value (see also Miiller & Biskamp (2000)). 

To our knowledge, no analytic model of comparable conceptual simplicity to those of 
Kolmogorov, Burgers, or Goldreich & Sridhar in their respective regimes presently exists 
for compressible MHD turbulence. In practice, elements of the latter two theories may both 
apply: at large scales, velocities are strong enough to drive shocks, while at small scales (away 
from shocks), velocity perturbations are small compared to the sound speed and magnetic 
perturbations are small compared to the mean magnetic field. However, the overall result 
cannot be a simple hybrid for the system as a whole, because there is no clean separation 
of scales when shocks produced by large scale motions transfer energy immediately to the 
dissipation scale. Results from lower-resolution numerical work indicate relatively steep 
power spectra (Mac Low & Ossenkopf 2000), while Cho & Lazarian (2002) have recently 
reported results showing anisotropy of the expected sense (see also earlier work in the weakly 
compressible limit by Matthaeus et al. (1996, 1998); Oughton et al. (1998)). The analysis of 
the present paper focuses on comparing turbulent spectra in strongly compressible models 
(A4 S = 5) with varying background magnetic field strength, based on simulations with the 
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highest resolutions (512 3 ) performed to date. 



3. Models and Parameters 

Our turbulent cloud models are obtained on a spatial grid of 256 3 zones (models A, 
B, C, and D) or 512 3 zones (models A2, C2). These models were studied in previous work 
(Stone et al. (1998), hereafter Paper II), which focused on the steady-state dissipation rate of 
turbulence. Models A-D here represent the same simulations with dissipation characteristics 
listed in Table 1 of Paper II, while models A2 and C2 are higher resolution counterparts. 

Each model begins as a cubic box of dimension L, initially filled with a stationary, 
isothermal plasma of uniform density po = p, and threaded by a uniform magnetic field B 
= B x. For each model, the initial magnetic field strength B is parameterized by (3 = 
c 2 /v\ = Cg / (B 2 1 '4tt p) , being proportional to the ratio of the thermal to magnetic pressure. 
We consider models representing a strong field case ((3 = 0.01; models A, A2), a moderate 
field case ((3 = 0.1; model B), a weak field case ((3 = 1; models C, C2), and a purely 
hydrodynamic case ((3 = oo, model D). 

We evolve each model in time using the ZEUS code (see Stone & Norman (1992a), Stone 
& Norman (1992b)) to integrate the compressible, ideal MHD equations 

| = -V-(PV), (1) 

at p 8np Anp 

— = Vx(vxB). (3) 

We adopt an isothermal equation of state P = c 2 s p, where the sound speed c s is constant 
in space and time. We view this as an acceptable approximation for modeling most of 
the gas within a GMC (see discussion in Ostriker et al. (1999)). We do not include any 
explicit viscous or resistive terms other than shock-capturing artificial viscosity in the MHD 
equations, and thus dissipation is purely numerical. 

Boundary conditions for all of our models are periodic. For a given model, (3 remains 
constant in time due to the periodic boundary conditions, which maintain volume averaged 
(5B)=0. The initial uniformity of B and p also implies a spatially and temporally constant 
mass-to-flux ratio, modulo the numerical reconnection effect that may occur if magnetic flux 
elements of opposite signs are advected into a single cell. 
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The initial velocity field of each model is v = 0, with turbulent driving implemented by 
introducing velocity perturbations 5v at time intervals At = 0.001t s (where t s = L/c s is the 
sound crossing time over the box). The perturbations follow a Gaussian random distribution 
with Fourier power spectrum 

(|5v(k)| 2 )ocA; 6 exp(-8A;/A; peafc ). (4) 

Peak driving occurs at k pea kL/2ir = 8 (X pea k = L/8) for all models except model Ck4, for 
which kp ea kL/27r = 4 so that it may later be used in a rescaling comparison with model C2. 
For all of our models, the energy input rate E due to driving, given by E/pqL 2 c^ = 1000, 
is constant. The driving is incompressive, with the constraint imposed that V • 5v=0. The 
perturbations are normalized such that no net momentum is added. This is accomplished 
by adding a spatially-uniform component of velocity to compensate for any net momentum 
associated with the sum of sinusoidal velocity increments, J p5v=0 . The resulting saturated- 
state energies associated with each model were presented in Paper II and will be further 
analyzed later in §3. 

Additional model parameters are provided in Table 1 (see also Paper II), including 
the turbulent Mach number M. s = <j v /c s and the Alfven Mach numbers A4a,o = &v/va,o, 
M.A,rms = &v/ v A,rms- Here the Alfven speed v A0 = B 2 /4np associated with the mean 
magnetic field is distinguished from v Arms = (B 2 ) /4irp, calculated via a mass- weighted 
average of the squared Alfven speed in each cell, |B + 5B| 2 /47rp, over the box. The velocity 
dispersion is calculated by a mass-weighted average over the box, a 2 = (|v| 2 p/p). The sonic 
Mach numbers are typically M. s ~ 5, while values for A4a,o vary from ~ 0.5 for the strong 
field models up to 5 for the weak field models. The total turbulent energies listed in Table 
1 are defined as E K = (1/2) f d 3 x\v\ 2 p and 5E B = (l/Sn) f d 3 x(\B\ 2 - _B 2 ), respectively. 
Dynamical time scales of interest include the flow crossing time tf(L) = L/a v and the 
Alfven wave crossing time tA(L) = L/va,o over the box, which are given by t s /M. s and 
tsfi 1 / 2 , respectively. Note that in Paper II, we frequently made use of the flow crossing time 
at the driving scale, tf(X pea k) = tf(L)/8 = t s /8Ai s . The values reported in Table 1 are for 
the saturated state at late times when the driving and dissipation have reached a balance. 2 

For the purpose of relating the dimensionless units of our simulations to physical units, 
we give the sound crossing time as t s = 53(L/10pc)(T/10K) -1 / 2 Myr, and the Alfven wave 
crossing time as t A = 7.6(L/10pc)(n^ 2 /10 2 cm- 3 ) 1 / 2 ( J B /10 / uG)- 1 Myr. The value of (3 for 



2 Valucs for the 256 3 models are averaged over times t — 0.25i s and 0.30i s . Due to computational expense, 
the 512 3 models were not evolved for as many time steps. Values for model A2 and C2 are given at times 
t = 0.10t s and t — 0.15i s respectively. These times are well beyond the time tf(\ pea k) ~ 0.025t s needed to 
allow for the development of shocks arising from bulk flow motions. 
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a given model determines the corresponding physical value of the initial magnetic field 
strength as B = 1.4 / uG / 3~ 1/2 (T/10K) 1/2 (n //2 /10 2 cm~ 3 ) 1/2 , while the sound speed is c s = 
O.lQkms-^T/lOK) 1 / 2 . The driving power input per volume is then given by E/L 3 = 
1.04 x 10~ 25 ergs cm _3 s _1 (n^ 2 /10 2 cm~ 3 )(L/10pc) _1 (c s /0.2kms -1 ) 3 . Because the current sim- 
ulations do not incorporate self-gravity, the length scale L is arbitrary. 



4. Results 



In this section, we present the results of our spectral analyses. §4.1 focuses on angle- 
averaged spectral properties; we evaluate inertial-range slopes for various components of the 
power, and assess how these depend on B . In §4.2, we turn to the issue of characterizing 
anisotropy in the power spectra. Spectral analysis was performed on the models using a 
Danielson-Lanczos fast Fourier transform (FFT) and computing the one-sided power spectral 
density (PSD). 

A question of some interest is the relative overall importance and scaling behavior of 
the compressive versus the shear modes of supersonic MHD turbulence. To compare these, 
we define separately the shear and the compressive component of the velocity PSD. These 
are computed in fc-space, after the FFT has been performed, as 



<w(k) = |kxv(k)|^ (5) 
^ 2 omp (k) = |k-v(k)| 2 (6) 

We investigate the behavior of these separate components in both §4.1 and 4.2. In 
Table 1, we give the volume-integrated energies associated with these two components, for 
all models. 



4.1. Angle-averaged power spectra 

To study power spectra as a function of length scale only (i.e. angle-averaged), we bin in 
spherical shells in k-space so as to generate PSD's as a function of k = k r = (A^ + A^ + A; 2 ) 1 / 2 . 
The resulting functions are the shell-averaged (average spectral energy per mode) power 
spectrum P(k) and can be related to the shell-integrated (total in a spherical shell) energy 
spectrum as dE(k) = Ank 2 dkP{k) /2. 

For model comparison, we will describe each of the resulting P(k) curves qualitatively 
in terms of four distinct ranges. These are, in order of increasing k (decreasing length): the 
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environmental range, driving range, inertial range, and dissipative range. The environmental 
range extends from the largest spatial scales (the computational box) to the scales at which 
the forcing is applied. For saturated, driven turbulence, the P(k) curve is nearly flat through 
these largest scales - i.e., there is near energy equipartition among modes. As k increases, 
there is a turnover shortly after k > k pea k, where k pea k defines the maximum in the forcing 
function (eq. 4), indicating the start of the inertial range. In the inertial range of k, the 
curve P(k) is often fit well by a power-law, with slopes varying roughly from -3.3 to -4.9, 
depending on model parameters. For the present simulations, the inertial range typically 
spans from kL/2n = 10 to 30 or 40, but is extended to 60 for the higher resolution models. 
In some cases, the choice of inertial range is somewhat subjective; the details of such cases 
are included in the discussion that follows. Finally, as k increases still, the curve turns 
steeper once again at what we identify as the (purely numerical) dissipation range. Through 
the dissipation range, slopes tend to range between -6.5 and -7. The dissipation range is 
determined by purely numerical effects, as the code does not explicitly include dissipation 
terms except for shock-capturing artificial viscosity in the MHD equations. This generic form 
for P(k) applies to both the kinetic and magnetic components of P(k), for both strongly 
and weakly magnetized models, as seen e.g. in Figures la,b and Figures 2a,b, respectively. 

The primary quantitative measure we evaluate from the P(k) in our models is the slope 
in the inertial range. Table 2 shows the exponent n of power-law fits, P(k) oc k~ n , in the 
inertial range for the PSD of the (specific) kinetic energy Pk(&) = v 2 (k)/2, the magnetic 
field energy Pb(/c) = 5B 2 (k)/87rp, the combination P tU rb(k) = -Fk(&) + Psik), the incom- 
pressive component of (specific) kinetic energy P s hear(k) = v shear(k)/2, and the compressive 
component of (specific) kinetic energy P CO mp{k) = v 2 omp (k)/2. The values presented in Table 
2 for the 256 3 models represent the analysis of an averaged data set obtained from two time 
snap-shots, t=0.25t s and 0.30t s (corresponding to 1.25tf(L) and 1.5tf(L) or 10tf(X pea k) and 
12tf(\ pea k))- The values reported for the 512 3 models A2 and C2 are at times t = 0.10£ s and 
t = 0.15t s , respectively. 

While many spectra from the models exhibit an easily identifiable power-law range, giv- 
ing typical errors in slope fits ±0.1, some important exceptions warrant detailed discussion. 
The PSD curves for density p (not shown) have no discernible inertial range for any of the 
models. Curves of Ps(k) typically have very little inertial range, and thus there is some 
ambiguity in choosing the boundaries of the range to fit; typically a range within kL /2n=10 
to 30 was used. The values cited were generated by averaging several fits around the center 
of the most likely power-law, resulting in slope fit errors of ±0.2. In several cases for P^{k) 
and Pshear(k), particularly those of the weak field model, there is a slight "break" in the 
middle of what we have identified as the inertial range, with shallower slope at larger k. In 
these cases, two approaches were taken to fitting the inertial range: (1) an single average 
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over the entire range and (2) fitting the two subranges on either side of the break separately. 
Results of the first method are cited in Table 2. Details of the second method are given in 
the discussion of the general results which follows. 

One striking effect is that the slopes characterizing the inertial ranges of Pturb(k), Pk(&), 
and P s hear{k) all appear to scale with (3, while those of the magnetic field and compressive 
velocity do not exhibit trends which are as definitive. Values for the PSD slopes of Pr(^) 
steepen as -4.0, -4.3, -4.7, -4.8 for 256 3 models of increasing (3 = 0.01, 0.1, 1.0, and oo, 
respectively. Steepening of the same sense is also evident for models A2, C2, with Pk(&) 
slopes -3.8 and -4.3, respectively. The presence of a magnetic field thus appears to reduce 
intermediate-scale dissipation of energy, in that the inertial range slopes become flatter as 
one goes from models of no or small magnetic field to those with greater mean field strength. 

Spectra Pb(/c) and P com p{k) are notably different from kinetic energy spectra. The 
slopes of Psik) vary only between -3.6 and -3.3 as one moves from the strong (A2) to weak 
(C2) field high-resolution models (with similar behavior at lower resolution). As the mean 
field weakens, the inertial range of Pb becomes progressively harder to identify, being almost 
entirely unrecognizable for the weak field models. The inertial range of P CO m P (k) is much 
broader than that of other variables, particularly for the strong field model where the power- 
law fits extend over nearly twice the range of the corresponding fits for P s / iear (/c). However, 
for moderate and weak- field low resolution models, the extent of the inertial range of P CO mp{k) 
is reduced again such that an average of fits is needed and the resulting errors in the indices 
reported are ±0.2. 

Notice that there are two main differences between the spectrum of low resolution mod- 
els and those for their high resolution counterparts. First, the indices in Table 1 for low 
resolution models tend to be roughly 0.2 steeper than those of the high resolution models. 
In some cases, this is due to the difference in identification of the inertial ranges. For the 
low resolution model, the inertial ranges are of lesser extent, which tends to result in fit 
being placed either closer to the turnover from the driving range or the turnover into the 
dissipation range; both cases lead to an artificial steepening of the fit. Second, several of 
the curves for P s hear{k) and -Ffc^), for both high and low resolution models, exhibit a slight 
"break" in the inertial range, as previously remarked. The consequence is that the average 
fit made through this break results in an apparent discrepancy between the indices reported 
for high and low resolution counterparts. 

The results of fitting either side of the slight spectral break within the inertial range 
with separate power-laws reveals that in fact there is agreement between the high and 
low resolution models, with some interesting caveats. The resulting indices of Pk(^) are 
4.0,4.3,(4.9:4.5), and 4.8 for models A, B, C, and D, respectively, with the notation (ni : n r ) 
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giving the indices of fits to the left and right of the break (when one exists) with individ- 
ual errors of ±0.1. Compare these to the indices found for high resolution counterparts A2 
(4.0:3.5) and C2 (4.9:3.9), again with individual errors of ±0.1. The fits to the left of the 
break within the inertial range are in agreement. The case of P s hear{k) is similar. Average 
and "left-right" indices n and (nj : n r ) of fits to P s hear(k) for low resolution models are 4.0, 
(4.5:4.0), (5.0:4.5), 4.9 for models A, B, C, and D, respectively. Again, comparing with the 
indices for the high resolution counterparts A2 (3.8:3.5) and C2 (5.0:4.0) reveals agreement 
to the left (smaller k) of the break within the inertial range. For Pb(/c), the two-part slopes 
are (3.5:3.7) for model A2, which may be compared with 3.6 for the "average" slope, and 
3.5 for model A. 

Note that when a break is apparent within the inertial range, it occurs well beyond 
the turnover that separates the driving range from the inertial range. Also, the portion 
of the spectra to the left of the break do not in any way mimic the input driving spectra. 
Although the current simulations do not have sufficient resolution to be definitive on this 
point, we speculate that the spectra to the right of the break (larger k) represent the true 
asymptotic regime, while the spectra to the left of the break (smaller k) represent a transition 
between driving and asymptotic regimes. The steeper slope in the transition range may 
reflect the envelope of wave distributions in the process of relaxing, via interactions, from 
the input spectrum to the asymptotic spectrum. This relaxation involves spreading in k- 
space accompanied by a decline in the peak amplitude. 

We now consider two examples to examine more specific effects which the magnetic field 
strength has on the form of the PSD curves. Figure 1 illustrates the strong magnetic field 
case, with the spherically binned PSD for a high-resolution 512 3 model overlayed with that 
of the lower resolution 256 3 model (models A2 and A respectively, both at time t = 0.10t s ), 
both of ,3=0.01. 

In Figure 1, there is clear indication of an inertial range between kL/2n = 10 and 60 
for model A2, showing the cascade of power from larger scale to smaller, dissipative scale 
at kL/2n > 60. The power- law fits to the inertial ranges of -Pr(^), P s hear(k), and Pb(/c) 
indicate scalings consistent with that of Kolmogorov turbulence (n = 11/3) or shallower, 
while that of P CO mp(k) is significantly steeper. Notice from Figure 1 that the lower resolution 
model (A) has a more limited inertial range than that of the high-resolution model. The 
doubling of the resolution in effect doubles the extent of the inertial range, because of the 
greater separation between the driving scales and the dissipative scales associated with the 
computational grid. 

This interpretation is supported by Figure 2, which illustrates the weak field model. 
Here we again plot P(k) curves for a high- resolution model, with k pea kL/2-n = 8. We overlay 
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these with the spectrum from a 256 3 model which was driven with k pea kL/2n = 4, rather than 
8 as before. Here, both models are shown at t — 0.15t s . For the 256 3 model, we have rescaled 
k by a factor of 2, and rescaled the magnitude by a factor of 10, in order to compensate for 
the change in k pea k and the difference in total energies of each model, respectively Note that 
the inertial ranges are of the same extent and that the curves practically line up, because 
the ratio of peak driving scale to dissipation scale are identical. 

While driving at lower wave numbers evidently provides a larger inertial range on the 
k > kp ea k side, our standard choice of k pea kL/2ir = 8 significantly larger than unity allows 
us to investigate inverse-cascade effects, that is, transfer of energy from the driving scale 
to the larger "environmental" scales. Evidence of the growth of power on large scales is 
presented in Figure 3, which shows time series of magnetic and kinetic power spectra from 
successive model snap-shots (times t = 0.05t s , 0.10t s , and 0.30£ s ). For both the strong- 
field and unmagnetized models, power at low k increases over time. The unmagnetized and 
magnetized models have similar, nearly flat, Pk(&) at low k, whereas Ps(k) rises slightly 
toward k = for the magnetized model. Thus, while indications of an inverse cascade are 
present, they may be principally due to compressibility. 

Note that Figure 3 further indicates that a steady-state saturation has been reached 
in the spectra for values of k > k pea k by time t = 0.05t s . Similarly in Paper II, Figure la 
shows that total energy saturation of forced models was also reached by time t = 0.05t s . 
The characteristic timescale for saturation of hypersonic (M. s > 1) turbulence is that which 
characterizes the bulk flow velocity at the driving scale X pea k, namely the flow crossing time 
tfi^peak) = tf(L)/8 = 0.025t s (L) for M. s = 5. This is because the relevant shock formation 
time scale is tf(X pea k), rather than the sound crossing time t s that would be needed for shock 
formation in a nonlinear wave steepening process. 

How important is the compressive component of the velocity field in the dynamics of 
turbulence? At moderate Mach numbers, it depends on the ambient magnetic field strength. 
In Figure 4, we show an overlay of P s hear(k) and P CO mp{k) for several models of differing field 
strength, illustrating the difference in magnitude and slope of the shear versus compress- 
ible components of the turbulent velocity spectra. Those shown in Figure 4 are the lower 
resolution 256 3 models, averaged over times t = 0.25t s and 0.30t s . 

Table 2 provides for several comparisons between models. For the strong field case 
(models A2, A), the compressive energy is an order of magnitude smaller than the shear 
energy, and has a steeper slope (-3.7 vs. -4.3 for the 512 3 models, or -4.0 vs. -4.6 for the 256 3 
models). As the magnetic field strength decreases, the difference in relative magnitudes and 
slopes of the shear and compressive power decreases until, for the pure hydrodynamic model, 
the difference in magnitude is about a factor of 3 and the slopes, -4.7 and -4.9, are almost 
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indistinguishable, given the errors in average fitting of P C om P (k) of model D. The magnitude 
difference in part simply reflects the fact that the forcing is incompressive. For weak field 
cases, the spectral slopes of the shear and compressive components are also very close (-4.4 
and -4.2 for the 512 3 model or -4.7 and -4.5 for the 256 3 model), given the afore mentioned 
details of fitting averages. 

The values we find for the total shear and compressive energies for each model, as 
listed in Table 1, indicate that E comp accounts for 24% of the total kinetic energy for the 
unmagnetized case (model D), 21% for model C, 19% for model B, and finally only 12% for 
the strong field case (model A). Comparison of models A, A2; C, C2 shows no dependence 
on numerical resolution of these ratios. 



4.2. Directional dependence of the power spectra 

We begin with a qualitative look at the turbulent anisotropy as seen in Figures 5, 6. 
These represent gray scale plots of (v/c s ) 2 and (SB) 2 / (Aitpc 2 ) from a 2D slice taken through 
the computational box (z=0) for the strong field case (512 3 model A2) at time t = 0.10t s . 
Notice that the structures are elongated along the direction of the mean magnetic field 
(horizontal axis). There are clearly stronger small scale variations in directions perpendicular 
to the mean field, suggesting more power in larger and in smaller ky. 

To study more quantitatively the global spectral anisotropy induced by the presence of 
the magnetic field, we have binned the PSD in annular shells, concentric about the axis of the 
ambient magnetic field direction x, so as to compare the PSD in the plane perpendicular to 
the field with that in the direction parallel to the field. The two-dimensional angle-averaged 
power spectra thus generated are functions of k±_ = (ky + k 2 ) 1 ^ 2 and k\\ = k x . The 2D contour 
maps of PSD (£;_|_, fey) provide a visual means by which to qualitatively characterize the effect 
of the magnetic field. 

We have computed two-dimensional power spectra for all our models, and find that 
they become more and more anisotropically distributed with increasing mean magnetic field 
strength, with more of the power concentrated along the directions perpendicular to that of 
the ambient magnetic field. 

Figure 7 shows contour plots of power spectra for the strong field case (f3 = 0.01, model 
A2, t = 0.10t s ). In addition to spectra of Pk(^_l, k\\) and Ps(k±, k\\), we also show the shear 
Pshear(k]_, k\\) and compressive P C omp(k±,k\\) velocity power spectra separately. The solid 
contours represent the levels of constant power, from log(P(fcj_,A;||)) = -3 (smallest k) to -9 
(largest k) in steps of 1. Three circular dashed curves at k=50, 100, and 150 are overlayed 
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to illustrate a perfectly isotropic relation between k± and k\\ for comparison. 

From Fig. 7, notice that the power remains concentrated at larger length scales along 
the magnetic field (small fcy) direction, while extending to smaller length scales in the two 
perpendicular directions (large k±). There is also a clear indication that the anisotropy is 
dependent on scale size in the P CO mp spectra (Fig. 7d). 

Compared to our results for the (3 = 0.01 model, we find that the power distributions are 
dramatically more isotropic with a decrease in field strength. Figure 8 shows contour plots 
for the power in the weak field case (f3 = 1.0, model C2, t = 0.15t s ). While the magnetic field 
and total velocity field are clearly isotropic, there is a small but real anisotropy along a 45 
degree angle between k± and k\\ for both the shear and compressive components of velocity. 
Note how the compressive field (Fig. 8d) has deficiency in power (extended contour) along 
the 45 degree line, while the shear field (Fig. 8c) has a surplus (contracted contour) along 
the same line; the total velocity obtained as the sum of these is, however, isotropic (Fig. 8a). 

To further investigate direction-dependent behavior of the power spectra, we have taken 
slices through our 2D contour plots along k± — 1 and k\\ = 1 to generate ID plots of P(k») 
and P(k±), respectively. For all our models, the qualitative morphological features of the 
directional P(k±) and P(k\\) spectra include driving, inertial, and dissipative ranges similar 
to the angle-averaged PSD curves already shown in Figures 1 and 2. 

For the strong-field model, there are clear differences in the directional power that reflect 
the spectral anisotropy evident in Figure 7; we illustrate the results in Figure 9. From Figure 
9, spectral anisotropy, as represented in the difference in magnitude of the overlayed P(k±) 
and P(k\\) curves, appears strongest in the cases of Pk and P s hear, with the scale dependence 
of this anisotropy clearly evident in the differing slopes characterizing the inertial ranges. 
The magnitude of the power with variations in the direction parallel to B is comparable to 
that of the power with variations perpendicular to the field only for the largest length scales. 
As k increases, P(k\\) falls off much faster than P(k±), such that the anisotropy increases 
with decreasing length scale. 

Within the driving range, the incompressible nature of the driving is clearly evident 
from the directional spectra in Figure 9, as seen in by the hump near k pea k which appears 
strongly in the P{k±) curve for v s h ea r, but not at all in the curves for v comp . Within the 
inertial range, both curves roughly follow a power law. It is difficult to identify a clear 
separation between the inertial and dissipative ranges for the P(k\\) curves. Note that from 
Figure 9d, for the largest scales there is isotropy in the distribution of compressive velocities, 
but as k increases, P CO mp{k±) becomes larger than P CO mp(k\\)- 

To construct quantitative measures of the anisotropy and to see how these vary with 
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changing magnetic field strength, we can compare indices of power-law fits through the 
inertial range of each of these P(k±) and P(k\\) curves, as done before with the angle- 
averaged power spectra. Table 3 shows the index n of the k~ n power-law fits for Pb, Pk, 
Pshear, and P CO mp for all models; in Figure 9, these fits are compared to the directional spectra 
for model A2. Note that such inertial range fits are only able to show an anisotropy in the 
strong field models. The previously mentioned subtle variations in the weak field models 
P comp and P sh ear (Fig. 8) are not seen here. 

There is a general trend that the P(k±) curves for the total and component velocity 
power spectra exhibit a steeper slope for the weak field than the moderate field case, and 
more so for the hydrodynamic case. This is especially seen for where the indices 

of the power-law fits are 3.5, 4.3, 4.7, and 4.8 for the (3 = 0.01, 0.1, 1.0, and oo models 
respectively. Note however that the reverse trend is evident for Pb (k\\), with spectral indices 
of 4.5, 4.0, 3.6 for the (3 = 0.01, 0.1, and 1.0 models respectively. 

The strong field in models A and A2 is responsible for the anisotropy in the power 
spectra as reflected in the relative slopes along k\\ and k± directions. The steepest profiles 
for any component of our models are for the strong field model along the direction parallel 
to the field, while the slope in the perpendicular direction for the same model is shallowest. 
For = 0.01 (model A2), P B (k\\) oc A;~ 4 ' 5 while P B {k±) oc A;' 3 ' . For the same model, 
Px{k\\) oc k^ 4 ' 7 and PK(k±) oc A; J 3 ' 4 . This is indicative of the fact that increasingly small 
amounts of power are transfered into the k\\ direction at large k. 

As we shall discuss further in §5, the exponents of the power-law fits to P(k±) and P(k\\) 
(for both Pk and Pb) can be related to a power-law scaling behavior between ku and k± 
that describes the spectral anisotropy (cf. GS). This can be quantified (see Cho & Vishniac 
(2000)) by asking, for a given power contour, what are the k± and ku intercepts? Starting 
with our 2D spectra, we compute a k»-k± relation by performing an interpolation of the P(k\\) 
curve into the P(k±) curve. The resulting function k\\(k±) provides information on how k± 
and fc|| are related across a given range of power contours. This scaling is not constant across 
all magnitudes of power, but is relatively constant within the inertial range. Figure 10 shows 
plots of k\\(k±) from the strong field model for P B , Pk, Pshear, and P comp - In the inertial 
range, the relation for Pb, Pk, and P shear is consistent with k\\ oc A;^ 3 or A; 3 / 4 ; the compressive 
velocity has quite different behavior, roughly k\\ oc k^ 8 . The scaling behavior for P comp holds 
over an extensive range, both at small k and for k within the expected dissipative range. 

We conclude this section with a discussion of one possible application of this work to 
observations. Namely, we propose a method for constraining the value of the mean magnetic 
field strength within observed cloud cores by analyzing the degree of anisotropy in centroid 
maps in molecular tracers. Consider a cloud core which has a nonzero component of the mean 
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magnetic field lying in the plane of the sky (xy-plane) and a relatively uniform density p (for 
the later to be true, the velocities must be sonic or below). Given spatially uniform molecular 
excitation, the centroid of an emission line at a position (x, y) in a projected map represents 
the mean value of the line-of-sight turbulent velocity, {v z )(x,y)= J v z (x,y, z)dz/ J dz. If one 
takes the Fourier transform of the 2D velocity centroid map, the result (v z )(k x , k y ) is identical 
to taking the k z = plane from the Fourier transform v z (k x , k y , k z ) of the full 3D data cube. 
In Figure 11a and lib, we compare the results of these two equivalent calculations, i.e. the 
FFT of the "centroid map" (v z )(x,y) and the k z = slice taken from the FFT of the full 
3D data cube of v z for the strong-field model A. The same anisotropy previously discussed, 
albeit with greater noise, is evident in projected maps. 3 

In Figure 11c and lid, we show that (v z )(k x , k y ) from the weaker-field models (B,C) has 
significantly lower anisotropy than does model A. Taken together, these results suggest that 
anisotropy in the Fourier transform of velocity centroid maps can provide a lower bound on 
the mean magnetic field strength in an observed cloud. From Table 1, A4a,o — (v 2 ) 1 ^ 2 / iu,o in 
model A is 0.5, while Ma,o m model B is 1.5. Here va,o is the Alfven speed in the mean field, 
which for projections along z is also the mean plane-of-the-sky field, va, p - Allowing for a 
factor of l/\/3 reduction for (v 2 ) 1 ^ 2 compared to (v 2 ) 1 ^ 2 , this implies that when (v 2 ) 1 / 2 /va, p is 
~ 0.3, anisotropy would be clearly evident (unless washed out by other effects such as strong 
but isotropic density variations), whereas when {v 2 ) 1 ^ 2 /va, p is ~ 0.9 or larger, negligible 
anisotropy is expected. 



5. Summary and discussion 

In this paper we present spectral analysis of driven, highly compressive (Ai s ~ 5) MHD 
turbulence from a collection of 3D numerical simulations whose parameters span a range 
believed to apply within GMCs. We investigate models that were the subject of previous 
work on dissipation of turbulence (Paper II), as well as new high resolution counterparts. 
Each of our models evolves an isothermal plasma in time subject to an applied incompressive 
turbulent driving of the velocity field. We parameterize the strength of the mean magnetic 
field B by (3 = c 2 /v\ and compare models of varying (3. Our main findings, and their 
relation to other work, are outlined as follows: 

1. The angle-averaged spectra for the total velocity and its shear and compressive 



3 We note, however, that anisotropy in velocity centroid maps can be washed out if the emissivity is highly 
nonuniform, as expected for any region having supersonic flow. Direct tests on our M = 5 models show, for 
example, that density-weighted velocity centroid maps have isotropic power spectra. 
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components (Pk, Pshear, Pcomp) all exhibit distinctive environmental, driving, inertial, and 
dissipative ranges. In some of the spectra, there is evidence within the inertial range (i.e. for 
wavenumbers not directly populated by forcing) of a slight break. In these cases, the portion 
to the left of the break may represent a transition region, and the portion to the right an 
asymptotic regime. For strong-field models, the spectrum of Pb also shows these distinct 
ranges. The spectra for Pb shows no inertial range for the weaker field models, while spectra 
of p show none for any model. 

2. Signs of a nonzero but weak inverse cascade toward modal energy equipartition 
(i.e. flat P(k)) is evident for both magnetized and unmagnetized models at k < k pea k- 
This develops within half a flow crossing time over the scale of the box (0.1t s , see Fig. 3), 
corresponding to four crossing times at the driving scale \ ea k- Observations of the turbulent 
spectra in clouds do not show such flattening towards large scales. Our work thus places 
serious constraints on the importance of turbulent driving at length scales f C i within 
GMCs, in comparison to inheritance of turbulence with an extended power spectrum from 
the process of cloud formation. 

3. For the strong field case (model A2), the values of the index n for the spectra of 
Pturb = (Pk + Pb), Pk, and P s hear are in the range 3.7 — 3.8, while the value for P B is slightly 
less steep, ~ 3.5. The velocity spectra averaged over the inertial range are thus consistent 
with the strong MHD turbulence model of GS, which has n = 11/3. In recent numerical 
work for incompressible MHD, Cho & Vishniac (2000) reported a ID spectrum consistent 
with n = 11/3, as did Miiller & Biskamp (2000) and Biskamp k Miiller (2000). 4 Maron 
& Goldreich (2001) report a shallower law consistent with n = 7/2 for their angle-averaged 
energy spectra; although they regard this as anomalous due to its inconsistency with GS's 
theory, we find it intriguing that our magnetic power spectra are similar. Interestingly, when 
we fit just the highest-A; portion of the inertial range, we also find n r = 3.5 for P K and Pshear 
for model A2. Cho & Lazarian (2002) also report a slope similar to 11/3 for Alfven modes 
in compressible MHD simulations, although their spectra are somewhat noisy because the 
inertial range falls at relatively low k in their models. 

4. For the weak field case (model C2), values n of the spectral index for P tur b, Pk, 
Pshear, and P comp are in the range 4.0 — 4.4, steeper than 11/3 and even exceeding the index 
4 associated with Burgers's spectrum. Although these measurements may be affected in 
part by spectral steepening toward the driving range, the slopes to the right of the "break" 



Both groups note that numerical schemes employing hyperviscosity show spectral flattening at large k 
from the so-called bottleneck effect; because ZEUS uses finite-differencing rather than a spectral method, 
our results do not show this bottleneck. 
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in the inertial range (e.g. 3.9 for P K ) remain steeper than those in model A2 (e.g. 3.5 for 
Pk)- Thus, a cloud permeated by a weak mean magnetic field - even if the total magnetic 
field strength is significant (cf. the total magnetic energy in Table 1) - would be expected 
to yield a distinctly steeper velocity spectrum than a cloud with a stronger mean field. Pg 
for the weak-field model has a much shallower nominal spectral slope, 3.3, but with large 
uncertainties due to fitting a limited inertial range. 

5. As the mean magnetic field strength decreases, the slopes of the spectra for Pturb, Pk, 
and P s hear all become steeper, implying that a stronger magnetic field inhibits dissipation 
somewhat at intermediate length scales. For a weak-field system, velocities are super- Alfvenic 
(leading to shock dissipation) over a larger range of scales than in a stronger-field system, so 
that the near-incompressive limit for which the GS model applies would be shifted to larger 
k. More generally, compressional effects become progressively stronger for decreasing field 
strength: the ratio of compressive energy E comp / E K =0.12. 0.19, 0.21, and 0.24 for models 
of P — 0.01, 0.1, 1.0, and oo, respectively. 

In part, the lower proportion of compressive energy in the strong-field models must 
reflect the tendency for energy injected into the Alfvenic cascade (which is purely shear) to 
remain there (cf. GS, Maron & Goldreich (2001), Cho & Lazarian (2002)). Because losses 
from compressive components are more rapid than from shear components, however, the 
steady-state ratio E comp /E to tai rnay differ from the relative losses through shocks vs. other 
dissipation (e.g. model A2 loses ~ 23% of the input power in shocks, while E comp /E total = 
0.07) 

6. Two-dimensional power spectra P(k\\,k±) in our strong-field model indicate global 
anisotropy, with spectral energy preferentially aligned along k± for Pb, Pk, Pshear, and P CO mp- 
Our weak-field models, however, do not show global anisotropy in their power spectra (but 
see below). These findings are consistent with the results of Matthaeus et al. (1998) and 
Oughton et al. (1998) that global spectral anisotropy (as measured using the mean anisotropy 
angle 6^) decreases as the ratio of fluctuating to mean magnetic field increases. 

7. Enhancement of anisotropy at large k is clearly evident in the two-dimensional power 
spectra P(k±,k\\), as well as in the cuts P(k±) and P(k\\), for our strong-field model. The 
global spectral anisotropy of Pturb, Pb, and P shear in strong- field models is consistent with 
the relation k\\ oc k 2 / 3 proposed by GS, while for P comp we find k\\ oc k 4 ^ 5 . 

The global spectral anisotropy scaling we identify in our strong-field model has previously 
been shown to describe local spectral anisotropy for simulated incompressible MHD turbu- 
lence by Cho & Vishniac (2000) and Maron & Goldreich (2001). As those works emphasized, 
the anisotropy predicted by GS is for k directions defined with respect to the local magnetic 
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field. If the large-scale field is strong in the sense that the turbulence is sub-Alfvenic (as for 
our models A, A2), then the local directions of field lines do not depart strongly from B ; for 
weak large-scale fields, on the other hand, local directions of B are uncorrelated with B (see 
Fig. 2 of Paper II). Thus, our finding of global spectral anisotropy only in our strong- field 
models is consistent with expectations, since only for models A, A2 does B correlate well 
with the local field. 

Discerning small-scale local anisotropy even when turbulence is super-Alfvenic on large 
scales is possible using second-order directional structure functions for coordinates defined 
with respect to the local magnetic field (Cho & Vishniac 2000; Maron & Goldreich 2001). 
Although we have not applied these specialized techniques here, our results for weak-field 
models are also consistent with the findings of Milano et al. (2001) that global spectral mea- 
sures may be isotropic even when local anisotropy is present as identified by other methods. 

8. The dependence of global velocity spectral anisotropy on the mean magnetic field 
strength suggests a potentially important observational diagnostic. Namely, observational 
evidence of turbulent anisotropy at a given spatial scale would imply that the value of the 
"mean field" Alfven Mach number - i.e. based on the plane-of-the-sky | (B) | for that spatial 
scale - is less than unity. In molecular cloud core regions where the density varies weakly in 
space, Fourier transforms of velocity centroid maps correspond directly to slices through the 
3D velocity power spectrum, so evidence of anisotropy is immediately accessible (see Fig. 
11). In regions where the flow is supersonic, convolution of the anisotropic velocity power 
spectrum with the isotropic density power spectrum tends to wash out anisotropy in velocity 
centroid maps, but more complex analysis techniques using the full position-velocity data 
cubes may still be able to ascertain the degree of anisotropy in the velocity power spectrum 
itself. Further investigation will show if such spectral diagnostics offer a practical indirect 
means of constraining the magnetic flux through a cloud core, and hence in determining 
whether it is magnetically sub- or super-critical with respect to gravitational collapse. 

9. Finally, we briefly comment on the relationship of our results to proposals presented 
in Boldyrev (2002) and Boldyrev et al. (2002) that the spectral slope for compressible MHD 
turbulence may be derived from a modified version of a She-Leveque analysis that pro- 
vides intermittency corrections to Kolmogorov scalings (see She & Leveque (1994); Dubrulle 
(1994)). Those works suggest that for a range of fractal dimensions for the dissipative struc- 
tures, the index of the power spectrum would lie in the range 3.74 — 3.89. For our fits over 
the whole inertial range, we find indices for Pk of 3.8 and 4.3, respectively, for models A2 
and C2, and indices 3.7 and 4.4 for just the shear component of velocity P s hear- Fitting just 
the high-A; portion of the inertial range for model A2 (C2), we find n r of 3.5 (3.9) for P K and 
3.5 (4.0) for P s hear- Thus, while our results for the strong-field model A2 are potentially con- 
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sistent with the Boldyrev et al proposal, the spectra of our weak-field model C2 are steeper. 
While detailed discussion of the reason for this difference is beyond the scope of this paper, 
we note that the basic assumption of a conservative cascade (implicit in the assumption that 
the slope of the third-order structure function £(3) = 1) is more applicable in the strong- field 
model that the weak-field model. It remains an open challenge to develop an analytic model 
that incorporates both direct shock dissipation and multi-scale energy cascades in a single 
framework. 

We are grateful to Charles Gammie for helping to initiate this project and to the ref- 
eree for a valuable and careful report. This work was supported in part by NASA grants 
NAG53840 and NAG59167. 
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Fig. I. — Power spectral density in strong field case ((3 = 0.01). Spherically binned P(k r ) of 
(a) P K (k) , (b) P B (k), (c) P sh ear(k), and (d) P comp (k). Both 256 3 (dotted curve) and 512 3 
(solid curve) resolution models are shown, with power-law fits (solid line) and indices (as 
labeled) for the 512 3 model in the inertial range. Heavy dashed line shows the input spectral 
driving shape (eq. 4) with arbitrary amplitude. 
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Fig. 2. — P(k r ) in weak field case (f3 = 1.0). Same quantities and fits as in Figure 1. 
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Fig. 3. — Time series of P(k r ) for (a) Ps{k) and (b) Pk(&) for the strong-field case (ft = 0.01, 
model A) and (c) Pk(&) for the unmagnetized case {(3 = oo, model D), at times equal to 
0.05t s (dotted curve), 0.1t s (dashed curve), 0.3t s (solid curve). 
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Fig. 4. — Compressive (dashed curves) and shear (solid curves) velocity power spectra for 
models (a) (3 = 0.01, (b) (3 = 0.1, (c) (3 = 1.0, (d) f3 — oo. Fits are made to the inertial 
range (solid line), as labeled. The corresponding total percentages of compressive energy 
(E cornp /E K ) are 12%, 19%, 21%, and 24% for (a)-(d). 



-27- 




Fig. 5. — Greyscale representation of (v/c s ) 2 in a slice (z=0) through the computational 
volume. Anisotropy is evident in the elongation of structures in the direction x of the mean 
magnetic field (left to right). To emphasize morphological features, logw (value) is shown 
from 0.0 (black) to 2.0 (white). 




Fig. 6. — Same as in Fig. 5, for (8 B) 2 / '(Att pc 2 s ) . 
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Fig. 7. — Strong field case (model A2, (3 = 0.01). Contour plots of P(k±, k») for (a) Pk, (b) 
Pb, (c) Pshear, and (d) Pcomp- Solid contours are log(P(k±,k\\))—-3 to -9, with the largest 
power being at the smallest k (bottom left of each pane). Dotted isotropic curves are included 
for comparison at k=50, 100, and 150. 
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Fig. 8. — Same quantities and levels as in Fig. 7, for weak-field case (model C2, f3 = 1.0). 
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Fig. 9. — Strong field case {(5 = 0.01, model A2). Overlays of P(fc||) (dashed line) and P(k±) 
(dotted line) for (a) Pk, (b) Pb, (c) P s hear, and (d) P CO mp- Fits are made to the inertial range 
(solid lines). 
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Fig. 10. — Measure of global spectral anisotropy for the strong field case {(3 = 0.01, model 
A2). Plots of fcy vs. k± for (a) P K , (b) Pb, (c) P s hear, and (d) P comp . Points of fey-intercept 
versus the /c^-intercept for a given power contour, as interpolated from the P(k±) and P(fc||) 
curves. 
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Fig. 11. — Spectral anisotropy in simulated observations as a constraint for magnetic field 
strength. A comparison of plots for (a) v z (k x , k y , k z = 0) and (b) (v z )(k x , k y ) for the strong 
field model A show that velocity centroid maps with uniform-density conditions evidence the 
same anisotropy as the underlying power spectrum. Much lower anisotropy in observations 
of (v z )(k x , k y ) are expected for conditions similar to the (c) moderate-field model B, and (d) 
the weak-field model C. 
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Table 1. Model Parameters 



Model 


P 


N 


Ms 


Ma,o 




Ek 


SE B 


-^shear 


Ecomp 


A 


0.01 


256 


5.11 


0.511 


0.477 


13.1 


7.35 


11.5 


1.55 


B 


0.1 


256 


4.86 


1.54 


0.986 


11.8 


7.14 


9.60 


2.20 


C 


1.0 


256 


5.08 


5.08 


1.67 


12.9 


4.13 


10.1 


2.77 


D 


oo 


256 


5.56 


N/A 


N/A 


15.5 


N/A 


11.8 


3.71 




1.0 


256 


6.61 


6.61 


1.75 


21.9 


6.42 


17.4 


4.76 


A2 


0.01 


512 


5.26 


0.526 


0.490 


13.8 


7.65 


12.4 


1.44 


C2 


1.0 


512 


5.16 


5.16 


1.55 


13.3 


5.05 


10.5 


2.81 



Note. — Energies in units pL 3 c^. 



-35 - 



Table 2. Comparative Spectral Slopes 



Model 


P 




P B (k) 




-Pshear (^) 


Pcomp (A;) 


A 


0.01 


3.9 


3.7 


4.0 


4.0 


4.6 


B 


0.10 


4.0 


3.6 


4.3 


4.3 


4.5 


C 


1.00 


4.3 


3.5 


4.7 


4.7 


4.5 


D 


oo 


4.8 


N/A 


4.8 


4.9 


4.7 


A2 


0.01 


3.7 


3.6 


3.8 


3.7 


4.3 


C2 


1.00 


4.0 


3.3 


4.3 


4.4 


4.2 



Note. — Exponents n for the power-law fits P(k) oc k~ n of energy component 
power spectra (PSD). Error estimates for the fits arc discussed in the main text. 
Note that the values of n corresponding to a Kolmogorov and Burgers spectra 
arc 11/3 and 4 respectively. 
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Table 3. Comparative Spectral Slopes of Directional Spectra 



Model 


P 


p B {k\\) 


P B (k±) 




Pa(k±) 


Pshear(k |) 


Pshear (&_L ) 






A 


0.01 


4.5 


2.7 


5.3 


3.5 


5.5 


3.5 


5.5 


4.1 


B 


0.10 


4.0 


3.6 


4.3 


4.3 


4.3 


4.3 


4.3 


4.3 


C 


1.00 


3.6 


3.6 


4.9 


4.7 


5.0 


4.8 


4.5 


4.5 


D 


oo 


N/A 


N/A 


4.8 


4.8 


4.9 


4.9 


4.7 


4.7 


A2 


0.01 


4.5 


3.0 


4.7 


3.4 


4.8 


3.3 


4.8 


4.0 


C2 


1.00 


3.3 


3.3 


4.3 


4.3 


4.4 


4.4 


4.2 


4.2 



Note. — Values of spectral indices n of power-law fits <x k "to P{k\\) and P(k±). 



